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Fast and slow rotators in the densest environments: a 
FLAMES/GIRAFFE IFS study of galaxies in Abell 1689 at z=0.183 
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ABSTRACT 

We present FLAMES/GIRAFFE integral field spectroscopy of 30 galaxies in the massive 
cluster Abell 1689 at z = 0.183. Conducting an analysis similar to that of ATLAS'^'^, we 
extend the baseline of the kinematic morphology-density relation by an order of magnitude in 
projected density and show that it is possible to use existing instruments to identify slow and 
fast rotators beyond the local Universe. We find 4.5 ± 1.0 slow rotators with a distribution in 
magnitude similar to those in the Virgo cluster The overall slow rotator fraction of our Abell 
1689 sample is 0.15 ± 0.03, the same as in Virgo using our selection criteria. This suggests 
that the fraction of slow rotators in a cluster is not strongly dependent on its density. However, 
within Abell 1689, we find that the fraction of slow rotators increases towards the centre, as 
was also found in the Virgo cluster 

Key words: galaxies; kinematics and dynamics, galaxies: clusters: individual: Abell 1689 



1 INTRODUCTION 

1.1 Galaxies and environment density 

Early type galaxies (ETGs), despite having masses and luminosi- 
ties that span several orders of magnitude, obey a number of tight 
phenomenological laws. These, collectively known a s "scaling re- 
lations", include the color-magnitude diagram (CMD:lBaumlll9 5 p ; 
IVisvanathan & Sandag3 Il977| - ^Sandage & Visvanathan"l978Jb'). 



the c olor- (J and Mg — a {Burstein et al. 1984; Bender et al. 
1993h relations and the fu ndamental plane dPressler et al.lll987l ; 



Diorgovskii & Davislll987l) . With their remarkably small scatter. 



no rotation, often exhibiting kinematically decoupled cores and 
misalignment between kinematics and photometry. The latter are 
flattened systems, compatible with rotational symmetry, where 
ordered, large scale rotation is important for the gravitational 
equilibrium. While overlapping with the existing dichotomy 
among ETGs, the new classification crucially crosses the boundary 
between Es and SOs, in that FRs populate both morphological 
classes. In deed ATLAS ^'^ (the volume limited follow-up s urvey to 
SAURON, ICappellari et al. 1 1201 Tal lEmsellem et al. found 
that many morphological Es are FRs. They suggest a new clas- 
sification paradigm bas ed on kinematics rather than morphology 



they impose strong constraints on the structure and evolution of 
ETGs. Morphologically, ETGs are either classified as ellipticals 
(Es) or lenticular (SO) galaxies. In the late 1980s, based on new 
and more accurate spectroscopy, Es were di vided into two grou ps: 
pressure supported and rotation supported (iBender et al l ll989h . It 
is thus particularly interesting to investigate how such dynamically 
distinct systems formed and evolved while still obeying the very 
same scaling relations. Environment certainly plays a major role in 

[ alaxy evolutio n, as witnessed by the morphology-density relation 
Dresslerlll980l) : systems in denser sun^oundings are more likely to 
be ETGs. 

The advent of integral field spectroscopy (IFS) has brought 
a wealth of information to the field. The SAURON survey dis- 
covered the existence of two kinematica lly distinct classes of 
ETGs, slow and fast r otators (SR and FR. lEmsellem et al]l2007l : 
ICappellari et al. I l2007h . The former are systems with little to 
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sitication paradigm bas ec 
jCappellari et al.l201 Ibl) . 



ATLAS '^'^ also presented the kinematic morphology-density 
relation (kT-E), linking the fraction of SRs (/sr) with the local 
number density of galaxies, /sr is insensitive to environment den- 
sity over 5 orders of magnitude, with a sharp increase observed 
only in the inner core of the Virgo cluster. ICappellari et al. I( l201ld) 
conclude by asking what would be measured in the denser environ- 
ments beyond the local Universe: does the fraction of slow rotators 
increase further or does it stay constant? 

Addressing this question would give further insight on the pro- 
cesses that drive galaxy formation and evolution, and is indeed the 
goal of this work. 



1.2 This study 

We exploited the unique capabilites o f the FLAMES/GIR AFFE 
multiplexed integral field spectrograph dPasguini et al.l20o3) at the 
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Very Large Telescope (VLT) to investigate internal kinematics of 
galaxies in the densest environment, so to extend the density base- 
line of the kT-E relation. After describing the observations in Sec- 
tion[2] we present the data reduction and analysis in Section[3] The 
results are presented in Section |4] followed by their discussion in 
Section|5]and a summary in Section|6] 



2 OBSERVATIONS 
2.1 Sample selection 

Abell 1689 is a mass ive galaxy cluster at redshift z — 0.183 
( IStruble&Roodlll999h . Its reg ular, concentric X-r ay contours sug- 
gest it is a relaxed system jLemze et al.l 120081) . An X-ray lu- 
minosity of Lx = 20.74 X 10**ergs~^ makes it consider- 
ably more lu minous than Coma, which has Lx ~ 7.21 x 
dEbeling e t al] Il996h and Virgo Lx = 8.3 x 
fSohringer et al.l 19941) . Assu ming 7-yr Wilkinson 



Table 1. A summary of the VLT/FLAMES GIRAFFE/IFU spectroscopy. 
The seeing was measured on site using telescope guide stars. 



10'*" ergs- 



10*'' ergs . 

Microwave Anisotropy Probe Cosmology jKomatsu et alj 1201 ll . 
f2m = 0.27, Q,A = 0.73, ho = 0.71) its comoving distance is 
741 Mpc, giving 1" per 3.0 kpc, so that GIRAFFE deployable in- 
tegral field units (see below) sample up to 1 Rc for most galaxies. 
GIRAFFE permits the observer to target 15 objects simultaneously 
and we chose to target 30 galaxies as a compromise between sam- 
ple size and integration time. Ou r selection was based on a cat- 
alogue from iHalkola et all ( |2006|) . and in order to gain the maxi- 
mum possible signal-to-noise ratio (SNR), we initially selected the 
30 ETGs with the highest surface brightness within Rc (includ- 
ing the brightest cluster galaxy). This sample was then subject to 
two practical constraints. We needed all of our targets to have high 
resolution HST imaging, which limited our choice to candidates 
in the innermost regions of the cluster. Physical constraints from 
the instrument (see Section |23t ruled out some targets in the most 
crowded regions, forcing us to re-select from a reserve list. This 
left us with 29 galaxies inside the HST field of view and one out- 
side (galaxy 20). 
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120 


0.60 
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120 


0.60 
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120 


0.60 
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27 May 2009 
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120 


0.60 
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25 May 2009 
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120 


0.50 
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26 May 2009 
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120 


0.60 


8 


2 


27 May 2009 
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120 


0.60 
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2 


27 May 2009 


4 


120 


0.80 


10 


2 


28 May 2009 


4 


120 


0.65 



a dedicated optical fibre. Alongside the 15 IFUs, each positioner 
plate also houses 15 sky fibres. These are fully deployable just like 
the former but carry only one lenslet. 

Since the magnetic buttons are larger (10") than the IFU field 
of view, they cannot be deployed closer than a minimum distance 
of 11" thus constraining the sample selection: galaxies closer than 
11" on the sky must be allocated on different plates, if at all. As a 
result, some targets lying in the most crowded regions of the clus- 
ter were omitted. We proceeded to divide the sample in two equal 
sets, with galaxies numbers 1 to 15 assigned to plate 1 and galaxies 
numbers 16 to 30 to plate 2. Each plate was exposed 5 times for 2 
hours, for a total of 10 hours exposure time per galaxy. 

We remark that, as detailed in Table [T] the seeing was compa- 
rable to the size of the lenslets (0.52"). This reduces the coiTelation 
between adjacent spaxels. 



2.2 Archival data. 

We used F625W band imaging from the HST Advanced Camera 
for Surveys, combined with g and / band GEMINI imaging. See 
lHoughtonet"aL|j2012h for a thorough description of these data. 



2.3 VLT data 

We present new data taken with the FLAMES/GIRAFFE spectrom- 
eter at the VLT Unit Telescope 2. The L612W filter gives a resolu- 
tion R of 11800 (the minimum allowed on the instrument) with a 
wavelength range of 5732-6515 A (4858-5521 A in the rest frame), 
which includes prominent absorption features of old stellar popu- 
lations (for comparison, SAURON has a wavelength range of 4800- 
5380 A). The observations were carried out between 24 May 2009 
and 29 May 2009, as detailed in Table [T] which also contains the 
observing conditions. The instrument provides 30 independent in- 
tegral field units (IFUs), deployable anywhere on the focal plane. 
These are arranged in two positioner plates, each hosting 15; each 
deployable IFU is positioned by a magnetic button, with an ac- 
curacy better than 0.08" and contains an array of 20 square mi- 
crolenses, each with a side of 0.52" on the sky. They are arranged 
in 4 rows of 6 (with 4 "dead" comers) for a total field of view of 
3" X 2". Each lenslet is then connected to the spectrometer with 
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3 DATA 

3.1 Data reduction 

We extracted the spectra using the standard ESO pipehn^E follow- 
ing the guidelines ESO olfero 

Each morning the telescope produces a number of calibration 
frames, including bias, lamp flats and arc lamp frames. To extract 
the spectra from the raw images we used the closest calibration 
available. The pipeline is organized into 9 "recipes", distinct appli- 
cations with a number of user configurable parameters: we used the 
default values unless otherwise stated. 

For each night we created a master bias out of the 5 raw frames 
provided. We used the method ksigma and the recipe gimasterbias, 
with the keywords ksigma.low and ksigma.high set to 3.0 to remove 
cosmic rays. 

We then proceeded to "fibre localization", (tracing the spec- 
tra on the chip). This is done using a set of 3 very high SNR lamp 
flat frames, in the recipe gimasterflat. At each spectral pixel on the 
frame the recipe determines the locations, in the cross-dispersion 
direction, of the light peaks corresponding to the centres of each 
fibre signal. A curve is fitted to each profile, and is stored as the 
trace shape. We used the standard unweighted summation to ex- 
tract the spectra (we set the keyword extraction.method to SUM). 
We set to PROFILE+CURVE the keyword biasremoval.method, as 
advised by ESO on the website, while the keywords fibres. spectra 
and fibres. nspectra were modified to take into account the oc- 
curence of both broken and unused fibres. The manual indicates 
that, using SUM, the contamination between neighboring spectra 
is less than 10% of the counts. The recipe also determines the 
pixel-to-pixel variation corrections and the fibre-to-fibre transmis- 
sion variations. The wavelength calibration was done separately for 
each night using the recipe giwavecaJibration. The resulting wave- 
length solution has an accuracy of 0.009 ± 0.033A and a resolution 
FWHM of 0.61 ± O.O7A. The science extraction was performed us- 
ing giscience. We set the parameters biasremoval.method to PRO- 
FILE+CURVE und flat.apply to TRUE. 



3.2 Data analysis 

3.2.1 Photometry 

We used g' and r' band Gemini imaging to create a catalogue of a ll 
galaxies in the observed region of the sky jHoughton et al.ll2012h . 
We applied cuts at / = 22 and in the related error (cr^/ , agi < 0.1 
mag). The resulting catalogue has been used to compute the num- 
ber density of galaxies (Section|4j3}, the cluster Luminosity Func- 
tion (LF) (Section [5TTJ and the cluster CMD (Sect ion |5TTt. HST 
imagi ng was used to determine de Vaucouleurs jde VaucouleursI 
1953h effective radii Rc (using the curve of growth method of 
Houghton et ^ l2012l) and ellipticities e, whenever this was pos- 
sible. In practice one galaxy (number 20 in Table |2j lies partially 
outside the ACS field of view, and takes its photometric parameters 
from the r' Gemini image. 

Since a large fraction of our sample is found in very dense 
regions, the surface photometry is often contaminated by that of a 
neighbour. Consequently the Rc va lues in Table[2]include a quality 
flag Q, as in lHoughton et alj j2012l). 

Following Cappellari et al. " 120071) we adopted the method 



^ http://www.eso.org/sci/software/pipelines/giraffe/giraf-pipe- 
recipes.html 

^ ftp://ftp.eso.Org/pub/dfs/pipelines/giraffe/giraf-manual-2.8.7.pdf 



of moments to determine ellipticities: after identifying the image 
isophotes we compute, for each of them, the position angle of the 
major axis PA, the ellipticity e and the surface area A. The elliptic- 
ity of the k-th isophote ek is defined by 



yi 



(1 - Efc 



(1) 



where Fi is the flux associated with the i-th pixel, and the coor- 
dinates {x,y) are drawn from the galaxy centre, with the x-axis 
along the photometric major axis. The sum is conducted on the set 
of all pixels comprised in the k-th isophote. We associate to each 
isophote an ellipse of area Ak equal to the isophote area, elliptic- 
ity ek and position angle PAk, and associate to it a radius defined 
by Rk = y^Ak/ir. SAURON and ATLAS'''^ based their classifi- 
cation on values computed at 1 _Ro- We therefore define Ec as the 
value of efe computed within the isophote of associated radius Ra. 
The results are listed in Table |2] We find them to be robust against 
changes in Ra, except for galaxy number 9 (Table |2j, which ex- 
hibits peculiar photometry, having an abrubt change in both e and 
position angle at a radius of ~ 0.5". 



3.2.2 Stellar kinematics 

Stellar kinematics were extracted using pPXF, a penalized max- 
imum likelihood algorithm developed by ICappellari & EmsellemI 
|2004 i). It fits the line-of-sight velocity distribution (LOSVD) in 
pixel space, by convolving a linear combination of stellar template 
spectr a with an LOSVD expressed by the truncated Gauss-Hermite 
series jvan der Marel & Franxlll993l;lGerhardlll993b : 



1 + ^ h^Hmiv) 



(2) 



where y = (i) ~ V)/a and the Hm. are Hermite polynomials. In 
practice however, our SNR was mostly lower than that (~ 60) re- 
quired to reliably measure the weights /13 and h4 so we fitted a 
Gaussian function, obtaining just V and a in the above expres- 
sion. While ATLAS^° team u sed MILES stellar template library 
Sanchez-Blazguez et al]|2006l) , its FWHM resolution of 2.54 A 



Beifiori et al.l 1201 11) was lower than that of our data (see Sec- 



tion [3T)> so we used the high resolution version (R = 400 00) of 
the ELODIE template library dPrugniel & SoubiranI I2OO1I) . with 
FWHM of 0.13 A. The two libraries span similar regions in the 
age-metallicity space; MILES reaches lower ages (~ 7 Gyr vs ~ 8 
Gyr) and includes some old, metal poor stars (Z ~ 1/100Zq), 
but these are not relevant when fitting ETGs, and the change of li- 
braiy is unlikely to introduce any significant bias when compared 
to ATLAS^'^ measurements. All ELODIE templates have a gap at 
A ~ 5414 A, so we cut the galaxy spectra at 5300 A. For each 
galaxy we computed a weighted average (with sigma clipping re- 
jection) of all the 20 spectra, and fed it to pPXF along with all the 
templates available in ELODIE. This resulted in « 15 templates 
being selected for each galaxy, and we use this subset to fit the in- 
dividual fibre spectra of the galaxy. We used formal errors derived 
by pPXF (we did not exploit the penalizing functionality of the al- 
gorithm). These are typically of the order of 15 kms~^ for V and 
17 kms"'^ for a, but they do not take into account the correlation 
introduced when log-rebinning. 
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Table 2. Our sample of 30 bright galaxies in Abell 1689. 



Galaxy 


Halkola 


RA 


DEC 


Mk 


Rc 







PA 


p(SR) 




loe 






(deg) 


(deg) 


(mag) 


(ai'csec) 


(1-3) 






(0-1) 




(Mpc~2; 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


17 


197.94862 


-1 


.15575 


-23.28 


0.65 ± 0.29 


2 


0.286 


20 


0.02 


0.412 ±0.092 


2.381 


2 


3 


197.94831 


-1 


.14250 


-23.16 


1.27 ±0.57 


1 


0.583 


163 


0.00 


0.585 ± 0.053 


2.408 


3 


12 


197.95400 


-1 


.13998 


-23.57 


0.71 ± 0.32 


2 


0.363 


47 


0.00 


0.611 ± 0.087 


2.279 


4 


34 


197.96433 


-1 


.14525 


-22.58 


1.15 ± 0.51 


2 


0.032 


235 


0.00 


0.455 ± 0.105 


2.333 


5 


22 


197.96239 


-1 


.15608 


-24.97 


2.84 ± 1.24 


2 


0.219 


113 


0.00 


0.252 ± 0.022 


2.928 


6 


47 


197.96764 


-1 


.16323 


-25.41 


2.39 ± 1.04 


2 


0.330 


266 


0.91 


0.116 ±0.022 


3.083 


7 


65 


197.97640 


-1 


.15380 


-22.99 


0.85 ± 0.38 


2 


0.108 


311 


0.02 


0.352 ± 0.100 


2.339 


8 


31 


197.95769 


-1 


.17172 


-24.86 


2.74 ± 1.20 


3 


0.090 


196 


0.40 


0.080 ± 0.026 


3.468 


9 


38 


197.95643 


-1 


.17477 


-24.28 


7.51 ± 3.28 


3 


0.077 


-41 


0.00 


0.428 ± 0.030 


3.748 


10 


43 


197.95238 


-1 


.18489 


-23.87 


1.69 ±0.74 


2 


0.049 


-157 


0.37 


0.158 ± 0.046 


2.658 


11 


37 


197.94754 


-1 


.18383 


-22.75 


0.74 ± 0.33 


2 


0.217 


-41 


0.00 


0.510 ±0.096 


2.489 


12 


32 


197.95437 


-1 


.17140 


-26.18 


13.92 ± 6.08 


3 


0.149 


-62 


0.95 


0.037 ± 0.024 


3.607 


13 


14 


197.94564 


-1 


.16432 


-23.64 


1.58 ±0.69 


2 


0.107 


44 


0.06 


0.214 ± 0.049 


2.696 


14 


6 


197.94071 


-1 


.16265 


-23.48 


0.62 ± 0.28 


2 


0.361 


88 


0.03 


0.444 ± 0.093 


2.582 


15 


2 


197.93615 


-1 


.15561 


-22.81 


0.93 ± 0.42 


2 


0.029 


83 


0.00 


0.356 ± 0.095 


2.130 


16 


7 


197.94359 


-1 


.15721 


-23.63 


1.81 ± 0.79 


2 


0.138 





0.00 


0.257 ± 0.043 


2.779 


17 


20 


197.95362 


-1 


.15161 


-23.04 


0.46 ± 0.21 


2 


0.558 


138 


0.00 


0.505 ± 0.092 


2.285 


18 


41 


197.96484 


-1 


.15373 


-23.55 


0.95 ± 0.42 


2 


0.174 


198 


0.00 


0.410 ± 0.094 


2.830 


19 


35 


197.96782 


-1 


.15578 


-25.13 


2.27 ± 0.99 


2 


0.246 


180 


0.76 


0.136 ±0.026 


3.321 


20 




197.98916 


-1 


.15270 


-24.38 


1.55 ± 0.19 


2 


0.040 


98 


0.00 


0.171 ± 0.064 


2.061 


21 


69 


197.97680 


-1 


.16486 


-22.56 


0.52 ± 0.23 


2 


0.278 


278 


0.00 


0.439 ±0.107 


2.504 


22 


75 


197.97635 


-1 


.18002 


-23.47 


0.86 ± 0.38 


2 


0.301 


193 


0.00 


0.626 ± 0.089 


2.622 


23 


61 


197.96103 


-1 


.17819 


-23.19 


0.57 ± 0.25 


2 


0.055 


292 


0.00 


0.416 ± 0.092 


2.822 


24 


70 


197.96526 


-1 


.19081 


-23.66 


1.19 ±0.52 


2 


0.245 


-124 


0.09 


0.290 ± 0.055 


2.411 


25 


60 


197.96108 


-1 


.18797 


-23.64 


0.85 ± 0.37 


2 


0.261 


-41 


0.00 


0.509 ± 0.092 


2.469 


26 


29 


197.95675 


-1 


.17549 


-24.91 


4.64 ± 2.02 


3 


0.116 


-26 


0.45 


0.120 ±0.027 


3.446 


27 


42 


197.95666 


-1 


.17153 


-24.94 


3.36 ± 1.47 


3 


0.154 


-53 


0.96 


0.107 ± 0.029 


3.660 


28 


8 


197.93336 


-1 


.18330 


-23.54 


0.83 ± 0.37 


2 


0.617 


-35 


0.00 


0.588 ± 0.094 


2.968 


29 


28 


197.95005 


-1 


.17060 


-23.63 


0.69 ±0.31 


2 


0.081 


96 


0.00 


0.445 ± 0.099 


3.060 


30 


4 


197.93953 


-1 


.16139 


-23.35 


0.93 ± 0.41 


2 


0.265 


70 


0.03 


0.370 ± 0.092 


2.594 



Note. — Column (1): galaxy ID number used throughout this work. Column (2): galaxy ID from lHalkola et all i200d) . Column (3): right ascension in 
degrees and decimal (J2000.0). Column (4): declination in degrees and decimal (J2000.0). Column (5): /f-band galaxy magnitude derived from the apparent r'- 
band magnitude and connected as detailed in Section|5] Column (6): Rc obtained with a curve of growth method and masking nearby objects, see Section l3.2.1l 
Column (7): quality of the Rc determination. A value of 1 is only given to the best fits. Values of 3 are assigned to objects with severe contamination. Column 
(8): ellipticity determined with the method of moments, inside the isophote of area tt R^. Column (9): position angle determined with the method of moments, 
inside the isophote of area vr R^. Column (10): probability that the galaxy is a SR, see Fig.|4] Column (11): Xr measured within the whole IFU field of view. 
Column (12): Mean surface density of galaxies inside the circle centred on the galaxy and containing its 3 closest neighbours. 



Due to the high spectral resolution and low SNR, we decided 
not to subtract the sky, but rather to fit it sim ultaneously with the 
stellar templates. Like IWeiimans et"aLl ( |2009|) . we provided pPXF 
with all the simultaneous sky spectra and let the maximum likeli- 
hood algorithm rescale them to best fit the data. 



4 RESULTS 

4.1 Kinematic maps 

The results of the kinematic extraction can be seen on Fig.[T] There 
are four frames for each galaxy, from left to right: high resolu- 
tion photometry (from either HST or GEMINI); low resolution GI- 
RAFFE spectrograph photometry; the extracted velocity map and 
the extracted velocity dispersion map. Above each galaxy we give 
the ID number; the celestial orientation is given by the black com- 
pass arrows (N and E). Comer spaxels and spaxels corresponding to 



broken/unused optic fibres are depicted in black. Although the spa- 
tial resolution is low, rotation can be clearly seen in some galaxies, 
while no such features are seen on others. 

We cannot detect kine matically decoupl e d core s (KDCs) and 
double a peaks (2-a) as in I Krai no vie et alj ilOl ih . because our 
spatial resolution is too coarse. If we try to detect SRs from the 
velocity maps by eye, we identify at most six: these are galaxies 
4, 8, 12, 20, 26 and 27. The overall fraction of SRs in the sample 
would then be 0.20, in l ine with what was found in the Virgo core 
dCappellari et al .1201 Ibl) . However we are subject to contamination 
from face on discs appearing as SRs, which increases /sr. 

We also highlight five more objects which, despite exhibit- 
ing large scale rotation, have misa ligned kinematic axes , a feature 
more common in SRs than in FRs jKrainovic et alj201ll ): these are 
galaxies 1, 3, 5, 9, 17 and 25. Galaxies 3 and 17 have very high 
ellipticities, and are thus unlikely to be SRs. Galaxy 5 has high ve- 
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Figure 1. Kinematic maps of tlie Abell 1689 sample. Each horizontal set of four images depicts one of the 30 galaxies in the sample. The first plot shows 
HST photometry (apart from target 20). Superimposed is the FLAMES/GIRAFFE footprint. The second plot is the reconstructed image from VLT integral 
spectroscopy, where each square is a spaxel, corresponding to a lenslet in the instrument. Superimposed is an isophote at either He, or the closest fraction that 
fits into the IFU footprint. The four black comers coiTespond to unused "dead" corners, while other black spaxels (seen in 11, 15 and 30) correspond to broken 
or unused fibres. The third and fourth plots depict the kinematic maps: velocity and velocity dispersion. The black compass arrows show North and East. The 
colorbar limits are given in kms~^. 
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locity dispersion, and also contains an inner disc (R = 1.5 kpc) in 
the HST imaging. 

4.2 Xr and kinematic classification. 

lEmsellem et al.l ( l2007h introduced the estimator A/j to mea- 
sure the projec t ed sp ecific angular momentum of galaxies and 
lEmsellem et al.l | |201 ll) further show how the combination of \r 
and ellipticity e conveniently captures the kinematic boundary be- 
tween SRs and FRs. \r is defined as 



i6I 



(3) 



where Fi, Ri, Vi and ai are the flux, distance from the galaxy cen- 
tre, velocity and velocity dispersion of the i-th spaxel; the sum is 
cond ucted over all spaxels inside some subset / of the IFU foot- 
print. [Emsellemet^ ([20o3, IIqIII) define Xr{Rc) as the value of 
\r computed inside the ellipse of area nRg (see Section D. 2. It . 
In our study however that ellipse may either not comprise enough 
spaxels to reliably measure \r{Rc), or be too large to fit inside 
the IFU footprint. Therefore we used existing SAURON data to esti- 
mate how our particular observing setup affects the measured value 
of \r. 



<l 



0.20p 
0.15 - 
0.10 - 
0.05 - 
0.00 - 
-0.05 - 
-0.10 - 
-0.15 - 
-0.20 z 



• 1 




* ! 
• 










; 6 27 




1 19 


• 



Re["] 

Figure 2. Simulated observation of SAURON data with 
FLAMES/GIRAFFE, at ^ = 0.183. AXr is the difference between 
\r(IFU) (computed from simulated obs ervations of SAURON d ata with 
GIRAFFE) and the value \r(Bc) given in lEmsellem et alj JIOOtI) . plotted 
agai nst R^. Slow/fas t rotato rs are denoted by red/blue dots, and classified 
as in lEmsellem et al.l )2007h (upper panel). In the lower panel we show the 
values of Re for our sample of galaxies in Abell 1689, where the symbol 
colour indicates SRs (red) and FRs (blue), according to a classification 
done using the \R-e plot (Fig. |4), as described in Section 14.2.21 Dashed 
vertical lines define regions over which we estimate biases and systematic 
errors (see Section 14.2.11 . 



4.2.1 Effect of pixelisation on \r. 

The original SAURON sample covers a wide range of ETGs types 
jde Zeeuw et alj[2002h , and its data is publicly availabl^. We use 
it to simulate observations with FLAMES/GIRAFFE, in order to 
determine how distance and reduced spatial resolution affect mea- 
surements of \r . For each galaxy we cr eated a kinematic model us- 
ing kinemetrjQ l lKrainovic et alj[200€l) : each model was then pro- 
jected to the distance of Abell 1689 and convolved with a seeing of 
0.8", before being "observed" with FLAMES/GIRAFFE. We cre- 
ated 10000 realizations of each model, adding Gaussian errors of 
15 km and 17 km for V and a respectively (Section r3.2.2| l. 

and proceeded to measure \r f or each of them. 

In their \r vs e diagram, lEmsellem et al. I ( l2007l I2OIII) plot 
values computed on the same SAURON spectrograph images, at the 
same spatial scale of 1 Rc. For the small galaxies in our sample 
however, R^ covers just a few pixels whereas the large galaxies 
have i?o larger than the field of view of the IFU. For this reason we 
cannot follow the ATLAS'^° prescription precisely. We therefore 
introduced Xr{IFU), defined as the value of A_r computed using 
all the available spaxels in the IFU field of view and show through 
simulation of the SAURON results that it is a satisfactory proxy for 
Xr{Rc). 

Fig. |2| shows AXr plotted against R^, where AA_r is de- 
fin ed as the difference bet ween Xr(IFU) and the value of Afl(iie) 
of lEmsellem etal J ( |2007|) . We can use this information to deter- 
mine the correction and the uncertainty that we need to apply to 
Xr{IFU) to obtain Xr{Rc). It is clear how our ability to recover 
the true value of Xr{Rc) improves with increasing R^. To make use 
of this information we separate the sample into three groups, based 
on Re (the divisions are at Re values of 1.15" and 1.70" which 

^ http://www.strw.leidenuniv.nl/sauron/ 

^ The IDL KINEMETRY routine can be found at 
|http://www.eso.org/~dkrajnov/idl/ 1 



naturally divide the Abell 1689 sample and are shown as vertical 
dashed lines in Fig.|2}. We find the following biases (mean offset) 
and systematic errors (dispersion): for galaxies with Re< 1.15" 
AXr = -0.06 ± 0.09; for galaxies with 1.15" ^Re< 1.7", 
AXr = -0.01 ±0.04, and for the remainder, AXr = 0.01 ±0.02. 

We corrected Xr{IFU) according to the biases measured, 
summing the systematic errors in quadrature to the random errors. 
This correction takes into account both the different spatial scale 
between Xr{IFU) and XR(Re) and the different spatial resolu- 
tion between Xr(IFU) and ec. In Fig. [3] we plot simulated values 
of Xr (IFU) against published values of Eo (from lEmsellem et al.l 
l2007h . Despite the aforementioned differences, there is little (10%) 
misclassification in our diagram, especially at high values of Re. 
We can calculate the probability distribution for the number of SRs 
(galaxies below the line defined by 0.31 x y/e and the green line 
in Fig. l3l lEmsellem et al.ll20lTI) . This is most easily done with a 
Monte Carlo approach. For each galaxy we assume Gaussian errors 
in Xr, truncated so that ^ Xr ^ 1 and sample 100000 times. The 
resulting probability distribution is Gaussian-like and we find 12.3 
±1.7 slow rotators, where the true value is 12. This justifies both 
our choice of Xr{IFU) to substitute for XR{Re), and the use of e 
computed at a different resolution and radius than Xr{IFU). 

4.2.2 Xr measurements and the statistical calculation of fsn 

In Fig. m we show the Xr{IFU) vs te plot for our Abell 1689 
data. Given the simulation in the previous section, the values of 
Xr{IFU) have been corrected by -0.06, -0.01 and 0.01 for galaxies 
with Re< 1.15", 1.15" s;i?o< 1.7" and Re^ 1.7" respectively. 
The errors include both the formal random error (from pPXF) and 
the systematic error (0.09, 0.04 and 0.02 for the three ranges of Re 
from the previous simulations). Given these errors we can calculate 
the probability distribution for the number of SRs, as done previ- 
ously for the simulated SAURON data. We find 4.5 ± 1.0 slow rota- 
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Figure 3. \b.{IFU) vs for our simulated FLAMES/GIRAFFE observa- 
tions of tlie SAURON sample of galaxies. The green lines, which separates 
SRs (below it) from FRs, has equation Ajj = 0.31 y/e. Red and blue dots 
denote SRs and FRs respective ly, according to the oiiginal SAURON classi- 
fication iEmsellem et al.l2007h . While \ji{IFU) has been measured on the 
redshifted and resa mpled data, the values o f e on the x axis are the original 
values published in lEmsellem et al J )2007l) . Despite the latter being mea- 
sured on much higher resolution than Xji(I FU), and at a different radius, 
the impact on the classification is low. Misclassified galaxies correspond 
either to red dots above the green line and blue dots below it. 



tors, corresponding to /sr=0.15 ± 0.03. Galaxy number 9, which 
has peculiar photometry and an uncertain value of e^, has no effect 
on the result, because its value of Xr{IFU) is greater than ~ 0.25 
(the maximum allowed for any SR) by more than 3a. 

lEmsellem et al. warn about using only Xr to assign a 

galaxy to either the slow or fast rotator class. The discrepancy be- 
tween the "by eye" classification and the classification here bolsters 
that warning. However, when studying galaxies beyond the local 
Universe, such a detailed analysis as was done by the ATLAS '^'^ 
team is unfeasable. We are thus forced to rely on a statistical ap- 
proach. 



4.3 Environment density 

For each galaxy i n the sample we compute d the local environment 
density following ICappellari et al. I( l2011bl) . We defined E3 as the 
number density inside the circular area centred around the target 
galaxy and encompassing three other galaxies. Density estimates 
were done using only valid targets in the catalogue of Section lSTI 



We applied a constant field correction of 0.49 gal arcmin^ , mea- 
sured averaging data from on e hundred 1 arcmin'^ fields from SDSS 
DR7 jAbazaiian et alj[2009t) . without correcting for cluster/groups 
contamination. For comparison, the minimum value found in our 
sample is 3.83 gal arcmin"^. In Fig.[5]we show /sr versus E3; we 
compar e it to the results of the ATLAS'"^ survey, and in particular 
to /sR l ICapDellari et al ]l2011bh . The Virgo core corresponds to the 
densest bin in ATLAS^^, with /sr = 0.25, double that typically 
found in less dense environments (/sr ~ 0.12). We probe environ- 
ments with values of logio E3 between 2.06 and 3 . 75 : the minimum 
is comparable to the core of Virgo, and the maximum is 1.7 dex 
higher. In this sense our work starts exactly where ATLAS fin- 
ished. We find a sharp increase in /sr with projected density, rang- 
ing from /sR= 0.01 in the least dense environment to fsK— 0.58 
in the densest environment. Errors due to misclassification, albeit 




Figure 4. X[i{IFU) vs tc for all target galaxies in Abell 1689. The green 
fine has equation = 0.31 \/e and separates fast rotators (blue dots 
above it) from slow rotators (red dots below it). Error bars are dominated by 
the systematic error (Section [4.2.21 . Notice that we corrected the measured 
value of \ii{IFU) by subtracting -0.06, -0.01 and 0.01, depending on 
for each galaxy (see again Section |4.2.2t . The solid magenta line represents 
the edge-on view of axisymmetric galaxies with f} = 0.70 X e, while the 
black dashed lines represent the trajectories of 6 of these galaxies (with 
e = 0.85, 0.75, 0.35) as their viewing angle goes from edge-on (on the 
magenta line) to face on (towards t he origin). For more in formation on how 
these models were constructed see lEmsellem et"ai] l l201lh 
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Figure 5. Fraction of slow rotators /sr over the ETGs population in 
A TLAS'*^, including Virgo cluster (green circles, soHd green fine), as given 
in lCappellari et al] )201 Ibl) . and fraction of slow rotators in our sample of 
Abell 1689 galaxies (red circles, solid red line). Numbers at the top are the 
total number of galaxies in that bin, with the same color code. The error 
bars for the Abell 1689 points represent the uncertainty in the slow rotators 
classification, as estimated in Section l4.2.2l The green square is the value of 
/sr that we measure resampling Virgo using our sample luminosity func- 
tion. The error bars are smaller than the marker size. The red square is the 
average fraction of slow rotators found in our sample. The lower, smaller 
circles and dashed fines are the fractions computed with respect to the total 
cluster population, for ATLAS'^'-' (green) and Abell 1689 (red); this is an 
estimate based on spirals and blue ellipticals counts. 
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large, show that the densest bin in Abell 1689 has a higher fraction 
of SRs than Virgo core (Fig.|5). The intermediate bin has a value of 
/sR compatible, within the errors, with both the field-group envi- 
ronments and overall Virgo cluster value but is however lower than 
the Virgo core, /sr in the least dense bin is lower than ATLAS'^° 
field and group values. 

However, considering the whole Abell 1689 sample, we find 
for an average value of logioEs = 2.77 that the SR fraction is 
0.15 ± 0.03 (red square in Fig.|5}, which is the same as the overall 
SR fraction in the Virgo cluster, when sampled in the same way 
(green square). Furthermore, both values are similar to the field 
and group samples in ATLAS'"^, suggesting little to no difference 
in /sR when it is averaged over the whole cluster. 



5 DISCUSSION 

5.1 Sample selection effects 

In order to assess the robustness of our result, it is important to 
study the relation between the sample of 30 galaxies and its parent 
population. Our sample selection, limited by both observing and 
instrument constraints, biases our study in different ways. In this 
section we discuss how these effects change /sr. 

Having in mind the ATLAS"'^ study of the Virgo cluster as 
a point of comparison, we determined the Abell 1689 K-band 
LF. We took th e r'-ban d catalogue (Section [3.2.1t and, following 
iHoughton et applied a k-correction to the GMOS r - 

band magnitudes. T he results have be en converted to Ks band (and 
Vega system), using Maraston] 1 2005 ) mod els, where we assumed 
an age of 10.4 Gyr 



Houghton et alj 2012h and passive evolution. 
We finally applied a cut at = —21.5 mag, thus matching 
ATLAS''^ parent sample selection. The result is shown in Fig. |6] 
(blue circl es), where we compar e it with the cluster RS (as deter- 
mined by iHoughton et"aLll2012l . red diamonds) and our sample 
(yellow squares). The Virgo ETG LF is also plotted (green trian- 
gles). 

Knowing the K-band magnitudes of our sample, we can show 
/sr as a function of magnitude. In Fig.[7] the value of /sr observed 
in our sample (red) is compared to the fraction of SRs over the 
ETGs population of Virgo (green). The two are, within the errors, 
remarkably similar; however we do not reach magnitudes beyond 
« —23 mag to probe the faint SRs. 



5.1.1 Red sequence bias 

Our sample falls entirely on the red sequence (RS), a property that 
was not sought after. We know that the RS does not necessarily 
trace the morphological ETGs population, as it can include red spi- 
rals and omit blue ellipticals. How many ETGs lying off the RS 
have we left out of our sample? A rich, relaxed cluster like Abell 
1689 comprises a very small fraction of spirals, particularly in the 
core. In fact the ratio between the RS LF and the cluster LF goes 
from 1 at the bright end to 0.70 at Mk= -22.5 mag. This means 
that including 10 "blue" galaxies in the faintest magnitude bins 
would remove the RS bias, leaving us with a color "fair" sample. 
These faint galaxies are more likely to be FRs (Fig.|7), so the bias 
introduced by selecting galaxies on the RS leads us to overesti- 
mate /sr. In fact, if we assume that these ten galaxies are all FRs, 
and that they are distributed spatially much like the observed tar- 
gets, we can determine the kinematic morphology-density relation 
for an unbiased sample (with respect to color) which we show in 
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Figure 6. Abell 1689 (blue circles) and sample (yellow squares) K-band 
LFs, obtained from GEMINI/GMOS g' band and r' band imaging, r' 
band magnitudes have been converted to Mk^ band as d escribed in Sec- 
tion[5] Red diamonds trace the cluster RS, as determined bv lHoughton et al] 
120121) . Green triangles represent Virgo ETGs LF, from ATLAS '^^ sureey 
(data available at www-astro.physics.ox.ac.uk/atlas3d/) 
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Figure 7. Fraction of SRs for Abell 1689 sample (red histogram) and in 
Virgo ETGs (green), as a function of magnitude. The number of SRs and 
the total number of galaxies in each magnitude bin are given at the top. For 
Abell 1689 the number of SRs has been rounded off. (Section [4.2.2t . 



Fig.|5](smaller red dots and red dashed line), with the relation from 
ATLAS^° (smaller green dots and green dashed line). The good 
overlap between ETGs and the RS in Abell 1689 causes /sr to 
stay the same, whether the fraction is computed against the RS or 
overall galaxy population. This is not true in a less relaxed, spiral 
rich cluster like Virgo (Fig.|5). 



5.1.2 Magnitude selection 

We know that in Virgo, /sr varies as a function of Mk (Fig.|7]l, and 
that the LFs of Virgo and Abell 1689 are different, in that Virgo is 
relatively richer in brighter objects (Fig.|6]. Since our sample is not 
fully representative of the Abell 1689 population, what bias does 
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this introduce in the measured value of /sr? A rigorous answer 
to this question is impossible, because we do not know if the SR 
LF varies as a function of redshift and/or environment. In particu- 
lar, Virgo is a small and dynamically young cluster, whereas Abell 
1689 is a massive, relaxed system. However, using a simulation, we 
can estimate what SR fraction we would measure in Virgo with the 
same selection effects present in our Abell 1689 sample. 

Let us assume that the SR fraction as a function of magni- 
tude is the same in Abell 1689 and Virgo (reasonable given Fig. [7} 
and that the galaxies were selected only on their magnitude and 
no other properties (not true, as discussed in Section 12.11 but rea- 
sonable given the substitutions required to comply with proximity 
constraints). Using the actual Abell 1689 LF, we drew random sub- 
samples from the ATLAS '^^ Virgo ETG population. These samples 
yielded /sR = 0.16±0.01,in agreement with the actual Virgo value 
of 0.16. Thus despite being biased towards brighter galaxies, we 
should measure the same /sr; this is because we sample down to, 
but not including, the faintest magnitude bin of ATLAS'^°, where 
/sR suddenly increases. 



5.1.3 Other factors 

We remark that the distribution of projected ellipticities of our 
sample is different from that observed in Virgo. Since the former 
is richer in round objects, and since SRs generally appear rounder, 
it follows that our sample could be biased towards higher values of 
/sr. It may well be that the clusters e distributions are different, in 
which case a higher fraction of round objects may increase /sr. 

Another possib le source of bias is the i ntrinsic shape of Abel l 
1689: according to lOguri et al.l (l2005h and ICorless et all l l2009l) . 
Abell 1689 is elongated along the line of sight, so that its mea- 
sured E3 is higher than what we would observe from another point 
of view. If the cluster length along the line of sight direction were 7 
times longer than the diameter of the sky projection, then the value 
of E3 observed would be ~ 7 times the unbiased value. Since the 
maximum reasonable value of 7 is ~ 3, in Fig.[5]log E3 is overes- 
timated by at most ~ 0.5, which does not significantly affect our 
results. 

Finally we remark that the corrections to Xr{IFU) that we 
derived in Section l4.2. 1 [ increase /sr ; had we applied no correction, 
we would have 3.8 ± 1.0 SRs, so an even lower value of /sr. 



5.2 General remarks 

Abell 1689 has a higher average density than Virgo, but the same 
value of /sr. In s ide the cluster, in agreement with the findings of 
ICappellari et al.l ( 1201 Ibl) . /sr rises with projected density. In the 
least dense region /sr is smaller than the ATLAS^° field/group 
value. Given the low number of galaxies per bin, we cannot rigor- 
ously claim that this is representative. However, a similar "deple- 
tion" i s observed in the outskirts of Virgo cluster (ICappellari et al.l 
l2011bl) . One explanation could be that massive SRs are driven by 
dynamical friction towards the centre of the cluster If these were 
originally distributed in the cluster like other galaxies, dynamical 
friction would reduce their orbital velocity and radius. Since this 
process is more effective on more massive galaxies, it would con- 
centrate SRs (more massive on average) with respect to other galax- 
ies. 



6 SUMMARY 

We demonstrated the use of FLAMES/GIRAFFE in IFU mode to 
perform a survey of 30 galaxies in Abell 1689 at z = 0.183. The 
data has sufficient quality and spatial resolution to classify the ma- 
jority of targets as either SRs or FRs. In summary: 

(i) we find, in agreement with ATLAS'^^ results, that SRs pop- 
ulate the high luminosity end of the LF; the SR LFs measured from 
the Virgo ATLAS^'^ sample and our Abell 1689 sample are identi- 
cal down to Mk ~ —23 mag. 

(ii) the fraction of slow rotators in our sample is /sr=0.15 ± 
0.03. If we apply the same selection criteria to all Virgo galaxies 
in ATLAS'*'^, we find the same fraction (assuming that the distri- 
bution of SRs with magnitude is the same in both clusters). This 
indicates that /sr is not affected by the average number density of 
the cluster. Both Abell 1689 and Virgo average /sr are in line with 
the ATLAS'^° value for field and group environments. 

(iii) the fraction of SRs increases towards the denser, central re- 
gion of the cluster. This is in agreement with what was found in 
Virgo, where SRs concentrate in the cluster core. This could be a 
consequence of dynamical friction, as SRs dominate the high mass 
end of the galaxy population. 

It is important to expand this study, both to further study Abell 
1689 down to lower luminosities and increase the number of ob- 
served clusters, to quantify the scatter in /sr . 
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